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ABSTRACT 

The energy spectrum of the cosmic microwave background (CMB) provides a powerful tool 
for constraining standard and non-standard physics in the primordial Universe. Previous stud¬ 
ies mainly highlight spectral distortions (ji-, y- and r-type) created by episodes of early energy 
release; however, several processes also lead to copious photon production, which requires a 
different treatment. Here, we carry out a first detailed study for the evolution of distortions 
caused by photon injection at different energies in the CMB bands. We provide detailed ana¬ 
lytical and numerical calculations illustrating the rich phenomenology of the associated dis¬ 
tortion signals. We show that photon injection at very high and very low frequencies creates 
distortions that are similar to those from pure energy release. In the //-era (z > 3 x 10^), 
a positive or negative chemical potential can be formed, depending on the balance between 
added photon energy and number. At lower redshifts (z < 3 x 10^), partial information about 
the photon injection process (i.e., injection time and energy) can still be recovered, with the 
distortion being found in a partially Comptonized state. We briefly discuss current and future 
constraints on scenarios with photon production. We also argue that more detailed calcula¬ 
tions for different scenarios with photon injection may be required to assess in which regimes 
these can be distinguished from pure energy release processes. 
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1 INTRODUCTION 

The immense potential of cosmic microwave background (CMB) 
spectral distortions as a unique probe of standard and non-standard 
physics in the early Universe has spurred renewed theoretical activ¬ 
ity (see Chluba & Sunyaev 2012; Sunyaev & Khatri 2013; Chluba 
2013a; Tashiro 2014, for overview). Commonly discussed spec¬ 
tral distortions are primarily due to energy release (Zeldovich & 
Sunyaev 1969; Sunyaev & Zeldovich 1970; Burigana et al. 1991; 
Hu & Silk 1993a), for example, because of decaying relic parti¬ 
cles (Sarkar & Cooper 1984; Kawasaki & Sato 1986; Hu & Silk 
1993b; Chluba & Jeong 2014) or the damping of primordial acous¬ 
tic waves at small scales (Hu et al. 1994a; Chluba et al. 2012b,a), 
which causes //-, y- and r-distortions. In the future, these signals 
may be detected with an experiment similar to PIXIE (Kogut et al. 
2011), a satellite concept which could provide a major leap forward 
in sensitivity over COBE/FIRAS (Fixsen et al. 1996) and may open 
a new unexplored window to the early Universe. 

While most theoretical studies focused on distortions caused 
by early energy release, spectral distortions can also be created by 
direct photon injection. Physically, the difference is that initially 
only the photon field is affected, whereas for pure energy release it 
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is assumed that the baryons and electrons are heated but no photons 
are directly added. For photon injection, part of the injected photon 
energy is converted to heat (through Compton scattering and pho¬ 
ton absorption), while the other part remains in the scattered initial 
photon spectrum. 

One inevitable distortion created by photon injection is caused 
by the cosmological recombination process through uncompen¬ 
sated atomic transitions in the hydrogen and helium plasma around 
z 10^ (see Sunyaev & Chluba 2009, for overview). At this late 
stage, Comptonization is already very inefficient so that the emit¬ 
ted photons just redshift due to the expansion of the Universe (e.g., 
Rubino-Martm et al. 2006; Chluba & Sunyaev 2006). The shape of 
the cosmological recombination spectrum thus directly reflects the 
dynamics of the recombination process. However, for photon in¬ 
jection at earlier stages (z <: lO'*), Compton scattering redistributes 
photons more strongly over energy, so that direct information of the 
initial photon injection process is (partially) erased. In addition, at 
low frequencies, double Compton (DC) and Bremsstrahlung (BR) 
emission/absorption become more efficient, so that detailed calcu¬ 
lations are required to interpret the signals. 

Other examples associated with photon injection/destruction 
at a wide range of energies are related to photon-axion conver¬ 
sion in the presence of primordial magnetic fields (Tashiro et al. 
2013; Ejlli & Dolgov 2014), graviton-photon conversion (Dol- 
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gov & Ejlli 2013), non-evaporating spinning black holes (Pani & 
Loeb 2013), superconducting strings (Ostriker c& Thompson 1987; 
Vilenkin 1988; Tashiro et al. 2012), gravitino dark matter (Lamon 
& Durrer 2006) or simply primary and secondary photons created 
in particle cascades from decaying or annihilating relics (Shull & 
van Steenberg 1985; Chen & Kamionkowski 2004; Slatyer et al. 
2009; Hutsi et al. 2009; Valdes et al. 2010; Slatyer 2015). 

For most previous estimates related to CMB spectral dis¬ 
tortions, the simplifying assumption that photon injection can be 
treated similar to energy release was applied; however, as we ex¬ 
plain here, this approximation is only justified in certain regimes of 
injection redshifts, Zi, and frequencies, v;. For instance, for injec¬ 
tion at very low energies, hv^ < lO^'^kT^, where Ty(z) is the tem¬ 
perature of the CMB blackbody at redshift z, photons are efficiently 
absorbed and converted into heat by the BR and DC processes be¬ 
fore Comptonizing. Similarly, at high frequencies, during the pre¬ 
recombination era (z > 10^), most of the photon energy simply 
leads to heating of the medium through electron recoil. However, 
at intermediated frequencies the problem is richer. 

Here, we carry out a first systematic study for the evolution of 
distortions created by photon injection'. We illustrate the rich phe¬ 
nomenology of the resulting distortion signal for single photon in¬ 
jection, providing both numerical results and analytical approxima¬ 
tions and estimates. In particular, we introduce the photon injection 
Green’s function (Sect. 2), which allows computing the distortions 
for more general photon injection scenario. We also briefly discuss 
constraints on single photon injection scenarios (Sect. 4) derived 
from COBE/FIRAS (Fixsen et al. 1996) and big bang nucleosyn¬ 
thesis (BBN) constraints (Simha & Steigman 2008). 

For the main calculations performed in this paper, we assume 
that photons are only injected in the CMB band, hv < 30kTy. In par¬ 
ticular, we restrict ourselves to energies below the pair-production 
threshold, above which significantly more processes become rele¬ 
vant (e.g., Svensson 1984; Zdziarski & Svensson 1989) also related 
to interactions with light elements (e.g., Kawasaki et al. 2005). 
In Sect. 3, we briefly discuss photon injection at higher energies, 
for which most of the energy is eventually converted into heat. 
However, since several soft secondary photons can be produced 
more detailed calculations seem to be required. For similar reasons, 
we highlight the necessity to perform more detailed calculations 
for different scenarios associated with photon injection to assess 
whether these can be distinguished from pure energy release. 


2 GREEN’S FUNCTION FOR PHOTON INJECTION 

For small spectral distortions, the photon Boltzmann equation can 
be linearized. The thermalization problem can thus be solved using 
a Green’s function method. The Green’s function for pure energy 
injection, Gth(v, z), was already discussed earlier (Chluba 2013b). 
The final spectral distortion caused by different energy release his¬ 
tories can thus be computed as (Chluba 2013b) 

r , d(Q/Pr) , 

Af(z = 0 ) = J G,,(y,z') dz'- (D 

Here d(QIPy)l dz' parametrizes the energy release history. The en¬ 
ergy release Green’s function is well approximated as a sum of 
p and y distortion with an additional temperature shift (Chluba 
2013b). For higher accuracy, the r-type (non-/j/non-y) distortion 

* Photon destruction can be treated as negative photon injection, and thus 
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(Chluba & Sunyaev 2012; Khatri & Sunyaev 2012a; Chluba 2013b) 
becomes important in the p-y transition era (1 O'* <zg3xl0^) and 
in principle allows distinguishing between different energy release 
scenarios (Chluba 2013a; Chluba & Jeong 2014). 

In contrast, here we discuss photon injection at different injec¬ 
tion frequency v; and redshift zi. Assuming that the injected pho¬ 
ton number is smail, we can again compute the final distortion us¬ 
ing a Green’s function approach, but this time we have to perform 
computations on a grid of injection frequencies and redshifts. This 
allows accelerating the computation of the resulting distortion for 
more general photon injection histories and also provides detailed 
insights that can be used to develop analytic approximations. 

In practice, we simply add photons to the photon distribution 
at Zi using a narrow Gaussian centered at Vi and then compute the 
final spectral distortion at z = 0 using CosmoTherm^ (Chluba & 
Sunyaev 2012). With this procedure we obtain the photon injec¬ 
tion Green’s function, Gin(v, v', z')> which describes the spectral re¬ 
sponse for different injection frequencies v' and redshifts z'. This 
Green’s function in principle has two parts, one that is similar to 
the energy release Green’s function and the other that arises directly 
from the scattered photon distribution; however, from the numeri¬ 
cal point of view it is easier to include both parts simultaneously. 
The spectral distortion can then be computed as 

A/v(z = 0) = ^ 

where 5 (v, z) is the photon production or source term. We define it 
such that at any moment the relative changes of the photon num¬ 
ber and energy densities are ANy(z)INy{z) = J[5(v, z)/fiv] dv and 
Apyiz)!Pyiz) = [UplkTyiz)] ^S(v,z)Av, respectively. Here, we set 
ttp = G”/G^' » 0.3702 with Gf = /- l)dx = k\f(k + 1). 
The Green’s function is thus normalized as 

^ jG,„(v,v',z')dv= ^a,p,(0). ^V,(0). (3) 

Below we discuss the dependence of Gin(v, v',z) on the injection 
frequency and redshift in the different distortion eras. We further¬ 
more give simple analytic approximations that can be used to com¬ 
press the information needed to describe the Green’s function for a 
range of scenarios, thereby accelerating the computation. 


2.1 Basic ingredients and assumptions 

2.1.1 Energy exchange by Compton scattering 

The energy exchange between electrons and photons is controlled 
by Compton scattering, which can be described using the Kompa- 
neets equation (Kompaneets 1956). The total Compton-y parameter 
for the evolution of the photon distribution between redshift z and 
today is given by 


P' kTy cr-iNeC 

Jo m^c'^ H(l-¥ z) 


dz»4.9x 10-"(l-^z)^ 


(4) 


where is the free electron number density and H(z) the Hubble 
parameter. The approximation is valid in the radiation-dominated 
era, assuming standard cosmological parameters. A similar param¬ 
eter, fe = / ^l^tTiN^cdt, can be defined by replacing the photon 
temperature with the electron temperature, however, for the stan¬ 
dard thermal history the difference only becomes noticeable at late 
times, z < 200 (see Fig. 1). 


^ www.Chluba.de/CosmoTherm 
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The dependence of yy on redshift is illustrated in Fig. 1 for 
the standard cosmology. The value of yy determines how efficiently 
photons are redistributed over energy. Comptonization is very effi¬ 
cient at z > fewx 10^, withy.^ exceeding unity by a large amount. In 
this regime, a yu-distortion (Sunyaev & Zeldovich 1970) is formed. 

After recombination (z < 10^), redistribution by electron scat¬ 
tering becomes negligible for photons with CMB energies. How¬ 
ever, the amount of Doppler broadening still reaches ^ 1% - 10% 
between z =: 10^ and == 10“*, an effect that for example is important 
for the calculation of the cosmological helium recombination lines 
(Rubino-Martm et al. 2008). At z S 10“*, a well-known y-distortion 
(Zeldovich & Sunyaev 1969), characterized by partial energy re¬ 
distribution of photons, is formed when the medium is heated. In 
both aforementioned cases, simple analytic approximation for the 
photon injection Green’s function can be found (see Sect. 2.2 and 
Sect. 2.3). During the yu-y transition era (10“* < z ;£ 3 x 10^), de¬ 
tailed numerical computations are required, as the photon spectmm 
can be found in partially Comptonized states (Sect. 2.4). 

We immediately mention that for very high-energy photons, 
well above the ionization thresholds of hydrogen and helium, the 
effective y-parameter does not drop abruptly after recombination 
occurred. This is because for high-energy photons it does not matter 
if the electrons are bound in atoms or free. Thus, if a high-energy 
photon interacts with a neutral atom during the post-recombination 
or reionization era, this leads to extra ionizations and production of 
non-thermal electrons and subsequent particle cascades for which 
more detailed computations are required (Chen & Kamionkowski 
2004; Slatyer et al. 2009; Hiitsi et al. 2009; Valdes et al. 2010). We 
shall leave a detailed discussion of this case to future work. 


2.1.2 Initial effective photon temperature 


At z = Zi, the CMB spectrum is assumed to be given by a pure 
blackbody at temperature T\. Immediately after the photon injec¬ 
tion, the number and energy densities of the photon field are 


Ny{z{) = <(r.) 


Pyizf) = p"(r,) 


1 -H 


^Ny 


ANy 


(5a) 

(5b) 


where xy = hvJkT\ is the dimensionless injection frequency. The 
injected number of photons is parametrized by ANyfNy. The num¬ 
ber and energy densities of the CMB blackbody are, respectively, 
given by AP‘(r) = (Snf c^)(kTIhf G” ^ 410cm-*(r/2.726K)3 and 
p^y'iT) = (SjThlc^)(kTlhYGf ^ 0.26eVcm-’(r/2.726K)\ so that 
pP' oc r“* and V” oc and p^* ^ a^'kTN^' ^ 2.70lkTN^\ We 
assumed injection of photons in a narrow line around Xy. After the 
injection, the comoving energy density of the photon-baryon sys¬ 
tem remains constant^ and due to the large excess in the number of 
photons over baryons, one can use Eq. (5b) to determine Ty, giving 


r, ^ 7’o(l -tz.) 


QfpXy ANy 


( 6 ) 


Here, we set the effective temperature of the photon field after the 
photon injection to Tyfz) = To(1-l-z), where To = (2.726±0.001)K 


^ This assumes that there is no other energy release process and that one 
can neglect the small energy extraction caused by the adiabatic cooling of 
baryons (Chluba 2005; Chluba & Sunyaev 2012) or by extra coupling of 
dark matter to baryons and photons (Ali-Haimoud et al. 2015). 



Figure 1. Dependence of the y-parameters, jy and y^, on redshift. After re¬ 
combination, the y-parameters drop strongly since the number of free elec¬ 
trons decreases exponentially. At late times, electrons drop out of equilib¬ 
rium with the photons so that ye < yy. Around zk = 5 X 10“*, we have 
.Ve — yy — O-l. The strength of line-broadening, Av/v ^ 2 jyy In 2, is also 
illustrated. After recombination, it becomes much smaller than Av/v ^ 10“^ 
and can be neglected for most applications. Photon drift towards lower fre¬ 
quencies through electron recoil is roughly given by Av/v ^ -{hv/kTy)yy. 


is the CMB temperature today (Fixsen et al. 1996; Fixsen 2009). 
Inserting Ty into Eq. (5a), we find 


Tj ^ Tod + zd 


1 -H 



3apXi \ 


1 ANy 


(7) 


which shows that for injection at xy = xq = (4/3)/ap x 3.6016 
we obtain Tj = Tofl + Zi)- In this case, the added number of pho¬ 
tons is exactly the amount required to increase the energy density 
of the photon field from T, to Tyizf) while having both the effective 
number and energy density of the photon field agree with a vanish¬ 
ing chemical potential. Explicitly this means that a full blackbody 
distribution can in principle be restored by simply redistributing 
photons over frequency (smearing out the narrow line). Following 
Hu (1995), we call this case balanced photon injection scenario. 
As we will see below, a second balanced photon injection regime 
can be found at very low frequencies, where photons are quickly 
absorbed and converted into pure heat (Sect. 2.2). 


2.1.3 Initial Compton equilibrium temperature 


In our computations, we simply use T^ = T\ for the initial electron 
temperature at zy. However, after the photon injection, T^ quickly 
approaches the Compton equilibrium temperature (Zel’dovich & 
Levich 1970) in the distorted photon field. Assuming d-function 
photon injection, we find 


J,eq ^ /fivVl +«)dv ^ ^ 


Ak f v^n dv 


1 -f 


im 


P An(x, z) 
4GP> 


r, 


, e ' i -I- 1 

1 + h:i -—7-4 


ttpXy ANy 


dx 


( 8 ) 


with Up!A X 0.0925, x = hvfkTi and the distortion of the photon 
occupation number. An, which for single photon injection initially 
reads An = Gj'x^^d/x - xf) ANy/Ny. This expression shows that 
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for Xj « 1 we have Tj'* « rj[l - {apX\l2)h.NyjNy\ < T^, while 
for Xj » 1 we obtain re‘' » 7’i[l + {apX^J4)h.NylNy\ > T,. For 
Xi X 3.83, we find « Ti, which at low redshifts (z < 5 x 10"^) 
defines the transition regime between net heating (x; > 3.83) and 
cooling (xi < 3.83) caused by photon injection (Sect. 2.3). 

Although by adding photons to the spectrum we increase the 
elfective temperature of the photon field, energy exchange with the 
electrons does not necessarily lead to net heating of the plasma. 
Photons injected at low frequencies up-scatter and thus cool the 
electrons (+baryons), while injection at very high frequencies leads 
to Compton heating through electron recoil. This energy exchange, 
respectively, manifests itself in negative and positive // and y distor¬ 
tion contributions due to subsequent interactions of the electrons 
with the CMB blackbody photons. Flere, we identify the regimes 
in redshift and frequency for which net heating and cooling of the 
thermal plasma occurs (see Fig. 10 for summary). 


2.2 Chemical potential era 


At high redshifts (z > 3 x 10^), Compton scattering is very efficient 
in bringing electrons and photons into kinetic equilibrium. In this 
regime, a pure /r-distortion (Sunyaev & Zeldovich 1970; Illarionov 
& Sunyaev 1975) is formed. Neglecting the effect of DC scattering 
and BR, the chemical potential and final CMB temperature can be 
determined from Eq. (5) using the conditions Ny^izf) ~ N^(Zi) and 
Py^(Zf) ~ Py(Zi) for a Bose-Einstein spectrum"* 


Afy (Zf) « <( 7 -.) 


pf (Zf) » pf (Ti) 


AT 

1+3 —— AI 2 Po 
1 + 4 ^ -M3P0 


(9a) 

(9b) 


where ~ 1.3684 and AI 3 » 1.1106. Introducing the constant, 
x"" = 4 AI 2 - 3 AI 3 » 2.1419, the solution of this system thus is 


Fo(zi) 

^ 3 

Apy 

A ANy 

f 


Py 

3”a7. 

= - -\^ 

L 



r 


.^Ny 


X 0.5185 X - 

■3.60ie 

> - 



L ‘ 


J Ny 

ATizO 

M\ 

2 ^Py 

Ml ANy Ml 

Ti 


Py 





1 ANy 



« 0.4561jUo + 

3”iv7 



[Xi - XoJ 


l ANy 


1 ANy 


( 10 a) 


(10b) 


where we used Apyjpy - Op X; ANyjNj and Xq = {4l3)/ap. These 
expression again show that for photon injection at x; < 3.6 a nega¬ 
tive chemical potential is expected (Hu 1995). In this case, photons 
up-scatter on average, cooling the electrons and baryons and caus¬ 
ing energy extraction from the CMB blackbody part to balance the 
total energy budget. For photon injection at Xj > 3.6, the opposite 
happens and a positive chemical potential is formed. 

We assumed that no extra energy in form of heat was released. 
Also injecting Apy/p^lj^ into the plasma (energy that heats the elec¬ 
trons or baryons directly without any change in the photon num¬ 
ber), we obtain an extra contribution (Sunyaev & Zeldovich 1970), 
Ph ~ 1.4006 Apy/py\^, to the resulting chemical potential. Thus, 
more generally, in Eq. (10) one can interpret Apyjpy as the sum 


Apylpy = Apylpy\^ + opX'^ANy/Ny. Theu Ae effective chemical po¬ 
tential and blackbody temperature, formed after a short time, are 


pS(zi) ~ 


3 Apy 
Py h 


3(yp 



ANy 


Ar(zi) _ M2.., , 1 

T, ~ 3 3 Ny ' 


( ll a) 

(llb) 


This shows that due to the extra heating the total effective chemical 
potential can in principle vanish or be positive, even if photon in¬ 
jection causes a net negative value for p. This also implies that for 
both photon injection and energy release, an interpretation of the 
final distortion constraints becomes more ambiguous. 


2.2.1 Including photon emission and absorption 

After the initial evolution of the photon distribution and electron 
temperature reaches a quasi-stationary state, the amplitude of p 
slowly reduces due to DC and BR emission. The final, observable 
chemical potential today is then expected to be roughly given by 

Po(z = 0)»pS(Zi)J'(Zi), (12) 

where JT(Zi) is the distortion visibility function (Chluba & Sunyaev 
2012 ), which can be approximated as JT(Zi) ~ ^ 

malization redshift Zp ~ 1.98 x 10® (e.g., Burigana et al. 1991; Hu 
& Silk 1993a). This factor takes the subsequent thermalization pro¬ 
cess for the p-distortion into account, leading to a strong suppres¬ 
sion of the resulting chemical potential caused by photon injection 
or energy release at Zi » Zp. More accurate expressions for the 
visibility function can be given (Khatri Sc Sunyaev 2012b; Chluba 
2014). Using the results of Chluba (2014), we find 

J'b*b(Zi) » 0.983 - Q.Q3U(zlZpf ^'>] (13) 

to represent the distortion visibility at 3 x 10® < Zi g 6 x 10® very 
well for the standard cosmology. We neglected relativistic temper¬ 
ature corrections due to Compton scattering, which become notice¬ 
able at Zi > 4 X 10® (Chluba 2014). 


2.2.2 Improvement at low frequencies 

While the above analysis captures the basic features for photon in¬ 
jection in the p-era, the situation is slightly more complicated at 
very low frequencies (x < lO^"* - 10“®). The most important dif¬ 
ference is that for photon injection at very low frequencies, DC 
and BR are very efficient and thus absorb extra photons before 
they can be transported towards higher frequencies by Compton 
scattering. In this case, photon injection again becomes similar to 
pure energy release, resulting in a small positive chemical poten¬ 
tial, although from Eq. (10) and (12) a negative chemical potential 
Po » -1.872Tbb(Zi) ANyINy is expected for Xj « 3.6. 

We illustrate this aspect in Fig. 2, where we present the 
numerical results for the final chemical potential obtained with 
CosmoTherm after photon injection at different frequencies and red¬ 
shifts. As expected, for xi 3.6, the net distortion vanishes, while 
p > 0 at Xi > 3.6. At Xi < 3.6, the chemical potential is nega¬ 
tive, however, instead of approaching po » -1.872Tb|j(Zi) ANy/Ny at 
Xi <K 3.6, p returns to zero and again becomes positive®. 


These expressions are found by integrating the Bose-Einstein occupation 

number, nsE = - 1]^ and expanding in NTIT\ and po. ® We confirmed this statement although it is not directly visible in Fig. 2. 


© 0000 RAS, MNRAS 000, 000-000 



Photon injection 5 



Figure 2. Final chemical potential, = 0), after single photon injection 
at Xj for different zj. For comparison, we show the simple approximation 
Mo(z = 0) » (3Q'p//t‘^) [x; - xq] J|J|j(Zi) AWy/A^y, with distortion visibility 
function J|^l^(Zi) ~ 0.983 Jl - 0.0381(z/zp)^ as well as the im¬ 

proved approximation, Eq. (15). 


To understand this aspect in more detail, let us consider the 
evolution of the injected photons in two steps. Compton scattering 
broadens the initial narrow line due to the Doppler effect. The av¬ 
erage photon energy is furthermore affected by Doppler boosting, 
electron recoil and simulated recoil, all processes that can be de¬ 
scribed using the Kompaneets equation (Kompaneets 1956). Since 
the number of photons does not change through Compton scatter¬ 
ing, the expressions given above are directly applicable. 

Once we include photon emission and absorption terms, a part 
of the injected photons can disappear before reaching the quasi¬ 
stationary ;u-evolution phase. Through the absorption process, their 
energy is immediately converted into heat and thus becomes equiv¬ 
alent to pure energy release. This effect can be incorporated by 
defining the photon survival probably, 'Ps(xi,Zi), which determines 
the fraction of injected photons, ANy/Ny\^ = 'PsiXi,Zi)tsNylNy, that 
sundve the evolution towards the initial /j-phase. Here, it is as¬ 
sumed that this transition occurs rather fast before the main ther- 
malization of the yu-distortion occurs. It is furthermore assumed that 
the DC emission caused by the injected photons is negligible, an 
approximation that we find to be valid in the discussion below (see 
Sect. 3). The number of injected photons that were absorbed and 
converted into heat then is ANy/Ny\^ = [1 - 'PsiXi,Zi)] tsNyINy. 
Since photon emission and absorption are most important at low 
frequencies, it is clear that 'Ps(xi,Zi) —» 0 for Xy —> 0. On the 
other hand, photons injected at high frequencies have a high sur¬ 
vival probability, so that PsiXi,Zi) —» 1 for xy —> oo. We will discuss 
the precise shape for Ps{x\, zi) below (Sect. 2.2.4). 

Even if the specific moment of absorption after the injection 
event determines how much energy the absorbed photon has, until 
that very moment a corresponding energy was already extracted 
from the thermal background (i.e., the bulk of CMB blackbody 
photons). Thus, the net energy injected by the absorption event, 
assuming that overall the evolution towards a yu-distortion occurs 
very rapidly, simply is ^pylpy\^ x apX\[l -Ps{Xi,z\)\tsNylNy. 
Hence, the total effective energy density that has to be thermal- 
ized is still tspyjpy = tspylpy\^ + tspylpy\^ = u^jX', ANyjNy, whete 
Apy/py\^ = apXiP^(X[,Zi)ANylNy. Howevet, the effective number 


of extra photons that need to be ingested is reduced to 


an; ANy 


Ps(Xi,Zi) 


ANy 


This means that Eq. (1 la) is modified to 


(14) 


t^oizd ^ 


3 tSPy 
N Py h 


3 q^/i 


^AN„ 


+ ^\x,- xoPAx„zd\^ 




J Ny 


(15) 


which shows that for Xj <sc 1 (or specifically xqP,,{x\,z\) < x\) 
photon injection again leads to a positive chemical potential con¬ 
tribution, while at high frequencies the expression in Eq. (11a) is 
recovered. Using the approximation Eq. (24) for'Ps(Xi,Zi), we find 
very good agreement with our numerical result for yuj (see Fig. 2). 
The approximation slightly degrades towards lower redshifts since 
the transition to the quasi-stationary yu-phase no longer is quasi- 
instantaneous, however, for estimates this expression suffices. More 
detailed estimates can directly rely on the numerical results ob¬ 
tained for different injection redshifts and frequencies. 


2.2.3 Photon injection Green’s function in the p-era 

Putting all this together, we have the simple approximation of the 
photon injection Green’s function in the yu-era 

G.„(y,y',z)^^[x'-xon(x',z)],7*(z)M(v)-t ^G(y) (16) 

where M(y) « G(y) [0.4561 - 1 /x] is the spectrum of a //-distortion 
with G(y) = (2/;yVc2) xe7(e^ - 1)^ and x = fiy/ytTo- We added a 
temperature shift term, oc G(y), since it is clear that the part of the 
initial energy that is not stored by the yu-distortion has been fully 
thermalized. We can determine A(v', z) using the normalization con¬ 
dition, Eq. (3), yielding 

T(y' ,z)~ Up (x' -\x’ -xaP^ (x',z)] JT*(z)) . (17) 

It is straightforward to confirm that for the single injection 
source function, 5(y, z) = hv ANy(z)INy{z)6(v - Vi)6(z - Zi), with 
(47r/c) f M(y)dy = (/rVSjp/To) and (4;r/c) / G(y)dy = 4py(ro) 
and Eq. (2) one obtains the correct chemical potential and temper¬ 
ature shift for the total spectrum. 


2.2.4 Photon survival probability in the p-era 


How can we determine the photon survival probability, which 
turned out to be so useful for the description of the distortion caused 
by photon injection? Clearly, this probability is determined by the 
competition between Compton scattering and DC and BR emis¬ 
sion. For this purpose, we can look at a simpler photon evolution 
problem, including only Compton scattering (at = Ty) and ab¬ 
sorption. The evolution of the distortion of the photon occupation 
number, An(x,y), is then given by (e.g., Chluba & Sunyaev 2012) 


dAn 



— An + An{\ + 2/2bb) 
ox 


Mx,y) . 

— - An, 

X n\j]y0y 


(18) 


where nbb = l/(e^ - 1), 9y - kTy/m^^c^ and A(x,>') determines the 
DC and BR emission rate. From Eq. (4), y(zi,z) = yy(Zi) - Vyiz) 
between two redshifts Zi and z < Zi- 

For photon injection at high frequencies (x > 1), the absorp¬ 
tion term can be neglected until the photon distribution is suffi¬ 
ciently smeared out over the whole frequency range through Comp¬ 
ton scattering. Once this regime is reached, the standard thermaliza- 
tion process starts, so that the survival probability until this phase 
is very close to unity (Fig. 3). 
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Figure 3. Survival probability for different injection frequencies and red- 
shifts in the fi-era. The curves were computed using CosmoTherm. The sim¬ 
ple approximation, f’s{x,z) ~ is also shown for z; = 2 X 10®. 


In the other extreme (x « 1 and y < 1), photons are so rapidly 
absorbed that Comptonization has not even started to become im¬ 
portant before most photons are already gone. In this regime, the 
survival probability is given by 


z) ' 
Z) • 


-Tabs(-^<Z) 


X 


' A(x,z') o-tAcC 


dz'. 


(19a) 

(19b) 


(4tr/3;r)e^Jdc, 


H{\ + z') 

Assuming that BR is negligible, one has A(x, z) 
where Jdc = f x‘*nbb(l + Ubb) dx x: 25.976 determines the DC emis- 
sivity of the CMB blackbody field (Lightman 1981; Thome 1981) 
and ax 1/137 is the fine-structure constant. With this, we find 


^abs,dc(-’ti z) ' 


1.2 X 10-2* 


-d+z)^ 


3.5 X 10-® 


yf- 


( 20 ) 


Strictly speaking, this expression is only valid at low redshift (z < 
5 X lO'*), when Comptonization is inefficient and BR absorption 
dominates (contradicting our assumptions), an aspect we will re¬ 
turn to below. In particular, it strongly overestimates the absorption 
optical depth during the fi-era, since Compton scattering was com¬ 
pletely neglected. 

At higher redshifts, Doppler broadening and boosting as well 
as CMB blackbody-induced stimulated recoil become important at 
low frequencies. Neglecting photon emission and absorption the 
evolution of a narrow line. A«(x, 0) = A 6(x - xif/x^ is given by 
(Chluba & Sunyaev 2008) 


An(x,y) 


-[\aUIXi)^yYlAy 




( 21 ) 


The maximum of the distribution thus evolves as xjfy) = x; e' 
(Chluba & Sunyaev 2008), which removes photons from the low- 
frequency part, where they have a large probability of being ab¬ 
sorbed. Because of stimulated scattering, the speed of this motion is 
smaller than what follows from the classical solution of Zeldovich 
& Sunyaev (1969), for which one finds x\(y) = x; e^^. Thus, a sim¬ 
ple correction to the absorption optical depth can be obtained using 


fabs(jC,z) : 


x 


A(xe>’<^-‘'*,z') e 

^2 H{\+z') 


dz'. 


( 22 ) 


which takes into account that the absorption probability drops as 


the photon drift towards higher frequencies, but does not account 
for the effect of Doppler broadening. For DC alone, this gives 

■rab.,dc(^. z) ~ -- Fa (V^)], (23) 

where Fyy(x) - (jy the Dawson integral. For « 1, 

one finds t*^^ j^(x,z) ~ Tabs,dc(^,z) oc however, for y^ > 1, 
one has T*^sdc(2^’^) * (2-2 x 10^^/x^) which increases signif¬ 
icantly slower than without the effect of up-scattering. Still, even 
this improved approximation strongly underestimates the survival 
probability at low frequencies, since line broadening caused by the 
Doppler effect, Av/v - 2 xjy In 2, also helps to transport photons 
towards high frequencies. 

It turns out that a simple approximation can be found using 
the quasi-stationary solution to Eq. (18). This yields 

n(^,z)*e-*®/x^ (24) 


where Xc(z) = xjA(x^.,z)IOy is the critical frequency at which 
Compton scattering takes over the evolution (Burigana et al. 1991; 
Hu & Silk 1993a). From Chluba (2014), we have 


4a- 


1 + z 


2x 10® 


1/2 


xf^ X 1.23 X 10-2 


1 + z 


2x 10® 


(25a) 

(25b) 


for the critical frequencies of DC and BR, respectively. To percent 
precision, the total critical frequency is x2 (xP*2)2 + (x®’^)^ (Hu & 
Silk 1993a), an approximation that is sufficient for our purposes. In 
Fig. 3, we show the comparison of our approximation for injection 
redshift Zi = 2 x 10®, finding good agreement with the numerical 
result. Overall, we find the approximation to be valid at the 10%- 
20% level for 3 x 10® < few x 10®. 

The solution, Eq. (24), also represents the frequency depen¬ 
dence of the chemical potential, fi(z, x) - /ro(z) (Sunyaev 

& Zeldovich 1970) during the quasi-stationary chemical potential 
evolution phase. This highlights that also there the factor 
should be interpreted as the probability to have a non-zero chemi¬ 
cal potential at a given frequency and redshift in the fz-era through 
the competition of emission and absorption and scattering. 

Equation (24) also allows us to estimate the frequency Xh <K 1 
at which the chemical potential is expected to become positive 
again. Erom the condition Xi « xo^s(-t, z), we find 


Xh(z) « Xc(z)/8.3 


(26) 


to reproduce the numerical results very well (see Eig. 10). Pho¬ 
ton injection around this frequency during the ^-era corresponds 
to a balanced injection scenario, with vanishing chemical potential. 
Below that frequency, photons are effectively converted into heat 
before reaching the high-frequency tail through Comptonization. 


2.3 Compton-y era 

We now discuss the effect of photon injection on the CMB spec¬ 
trum during the y-distortion era (z < (O'*). At these redshifts, the 
y-parameter is already rather small, y^ ;£ 10-2 (see Eig. 1), so that 
any smearing of photons over frequency or direct energy exchange 
with electrons is quite limited. At low frequencies, the photon emis¬ 
sion and absorption process is furthermore dominated by BR. 
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2.3.1 Photon injection after recombination 


As Fig. 1 illustrates, for scenarios with late photon production, af¬ 
ter recombination finished (z < 10^), one can practically neglect 
Compton scattering, unless the initial photon energy is very large 
so that electron recoil, Av/v =: —xyy, or ionizations of neutral atoms 
become significant. Compton scattering remains negligible even if 
one accounts for the reionization and structure formation processes 
(e.g., Hu et al. 1994b; Cen & Ostriker 1999; Refregier & et al. 
2000; Oh et al. 2003) or heating due to magnetic fields (Jedamzik 
et al. 2000; Sethi & Subramanian 2005; Kunze & Komatsu 2014; 
Chluba et al. 2015), the former of which could on average increase 
the effective y-parameter to the level of 2 x 10“* from small haloes 
(Hill etal. 2015). 

Here, we do not consider photon injection at energies above 
the hydrogen ionization threshold, Vh - 5.8 x lO^'/fl -l- z), around 

the Lyman-a excitation frequency, Xh, ex - 4.3x 10 "^/(l-l-z), or at the 

corresponding energies of helium. This would affect the ionization 
history (Peebles et al. 2000; Chen & Kamionkowski 2004; Pad- 
manabhan & Finkbeiner 2005) and produce CMB spectral distor¬ 
tions from the reprocessing of photons and heat by atomic species 
(Chluba & Sunyaev 2009; Chluba 2010). In this case, a more de¬ 
tailed treatment is required (Slatyer et al. 2009; Hiitsi et al. 2009; 
Valdes et al. 2010; Slatyer 2015); these considerations are also rel¬ 
evant when deriving CMB anisotropy constraints on annihilating 
dark matter particles (e.g., Galli et al. 2009; Hiitsi et al. 2009, 2011; 
Planck Collaboration et al. 2015). Similarly, we neglect other inter¬ 
actions with atomic and molecular species at late times (e.g., 21cm 
absorption, Lyman-Werner bands, etc.). 

In this case, the photon distribution only evolves through BR 
emission and absorption. The kinetic equation for the photon occu¬ 
pation number then reads 


dn 

dr 


Abr e 


[1 - « (e'' - 1)] -I- x). 


(27) 


where = hvjkT^ and Abr determines the BR emissivity. We 
also included a possible photon source term, j(t, x), and defined 
the Thomson optical depth r = f crjN^cdt. As this equation shows, 
even for the CMB alone a distortion arises if the electron temper¬ 
ature differs from the CMB temperature, Ty, when free-free 
emission tends to bring the CMB spectrum into equilibrium with 
the electrons at low frequencies, enforcing n 1/Xe at x^ <K 1. 
This is independent of any extra photon injection and can be stud¬ 
ied separately (Hu 1995; Burigana et al. 2004; Chluba & Sunyaev 
2012; Trombetti & Burigana 2014), so that we neglect it below. 

The evolution equation for a distortion to the CMB blackbody 
then becomes 


dAn 

dr 


ABR(T,Xe)(l-e-G 


An + s ( t , x). 


(28) 


Between z = 0 and Zi, this equation has the simple solution (com¬ 
pare also, Hu 1995) 


An(x,z = 0) ~ An(x, Zi)e 


Tg(x,z) = 


r 


■f' 


tb ( x , z ') Jf) dz 


ABR(z,Xe)(l -e g-jN^cdz 

xl H(l+z)’ 


(29a) 

(29b) 


where we introduced the source function, s(z, x), with respect to 
redshift. At very low frequencies, BR absorption efficiently de¬ 
stroys photons, converting their energy into heat, which then reap¬ 
pears as a small y-distortion at high frequencies. However, since 
this occurs only at very low frequencies (x <K 10“'*), a lot of low- 
frequency photons need to be injected to have an appreciable effect. 



Figure 4. Free-free absorption optical depth at different redshifts and fre¬ 
quencies X. At low frequencies, the scaling is Tff ^ f (z) ln(2.25/x)/x^. For 
X ^ 10“'* (= 6 MHz), the Universe becomes transparent (rff ^1) around re¬ 
combination. For X ^ 10“^ and x ^ 0.01, this transition happens at z ^ 1700 
and ^10^, respectively. 


The solution, Eq. (29), shows that the Green’s function in the 
considered regime can be expressed as 


Gin(v,v',z) : 


- v') 
4;r 




X a„. 


(30) 


where Y{v) x G(y) [x coth(x/2) - 4] is the y-distortion. The last 
term accounts for the aforementioned small heating of the medium 
by BR absorption, which usually is negligible. At z < 200, when 
matter and radiation thermally decouple, the last term is further 
suppressed. We furthermore mention that for very small, but non¬ 
zero, y-parameter one can alternatively replace the d-function using 
the Compton scattering kernel (Sazonov & Sunyaev 2000). 

The free-free absorption optical depth can be calculated using 
the results from CosmoRec (Chluba & Thomas 2011) for the ion¬ 
ization history and approximations for the free-free Gaunt factors 
from Itoh et al. (2000). Alternatively, at low temperatures, one can 
approximate the free-free Gaunt factor as (Draine 2011) 


gff(jCe, Te) ~ 1 + In 


5.158 + exp 


a 

n 


In 


2.25 



(31) 


where Z, is the charge of the nucleus. For even simpler estimates, 
one can use gff(Xe) ~ ln(2.25/Xe) at Xe < 0.37 (Hu 1995). Ne¬ 

glecting helium, with Abr = {aAll27t'j^)Np6j'’^^ gfiix^), where 
Te ~ 2.426 X 10“*** cm is the electron’s Compton wavelength. 
Alp Ve the free proton number density and 6^ = kTdm^c^ x 9y, 
we finally obtain Abr ~ 1.9 x 10“’Ze(l -l- z)^’^^ ln(2.25/x) at 
X < 0.37, where is the free electron fraction. The free-free op¬ 
tical depth is thus roughly given by, Tff(x, z) ~ f (z) ln(2.25/x)/x^, 
where F{z) is a single redshift-dependent function. 

We illustrate the numerical results for Tff(x, z) in Fig. 4 for 
several values of x. In the post recombination era (z < 10^), sig¬ 
nals produced at x > 10““* (= 6 MHz) are not significantly at¬ 
tenuated by free-free absorption. For percent-level precision, one 
does, however, need to include the effect of free-free absorption at 
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Figure 5. Photon injection Green’s function for injection at low redshifts, zi 5 x 10“^. In purple we show the analytic approximations for the Green’s function 
using Eq. (35) for jc < 1 and Eq. (39) for higher frequencies. Notice that for several cases the approximation fully covers the numerical result. Usually, the 
high-frequency y-type distortion contribution is relatively small, unless photons are injected at xi ^ 1 /yy, for which heating through recoil becomes important 
(e.g., see panel for xi = 15). To convert between x and observing frequency today we may use x ^ 0.017 (v/GHz). 
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Photon injection 9 


X < 0.1 (= 6GHz) at z 10^ - 10^, an aspect that, e.g., is im¬ 
portant for the low-frequency hydrogen and helium recombination 
spectrum (Chluba et al. 2007a). 


2.3.2 Photon injection at £ z ;£ 10'* and x < I 

At redshifts 10^ < z £ 10*, the total y-parameter can reach the 
percent level (see Fig. 1). In this case, line broadening through the 
Doppler effect can be as large as == 10%, but no signihcant net en¬ 
ergy exchange between the injected photon distribution and elec¬ 
trons occurs for x, « l/yy. Thus, the Green’s function for this 
regime has two parts, one that is sourced by the absorption of pho¬ 
tons at low frequencies, where BR is efficient and causes a small 
y-distortion, and the other part related to the slightly scattered and 
broadened injected photon distribution plus a smaller y-distortion 
due to energy exchange. Both of these aspects can be approximately 
treated independently. 

At low frequencies, BR absorption effectively destroys pho¬ 
tons, and the photon survival probability is given by 

n(x,z)~e-^f'«>, (32) 


with Tff(x,z) ~ F{z)\n{2.25lx) x^^ from Eq. (29b). It is straight¬ 
forward to determine the frequency at which most (= 99%) of the 
injected photon energy is absorbed and converted to a y-distortion. 
At 10* < z :£ 10*, we find this for x - fewx 10“* in agreement with 
our detailed computations. 

At slightly higher frequencies (0.01 < x < 1), we can use 
the solution Eq. (21) to account for the effects of electron scatter¬ 
ing (Doppler broadening, Doppler boosting and stimulated scatter¬ 
ings). In this regime, the average energy of the photon distribution 
increases like Apy(yj)lpy = (Apylpy)e^^T. The energy required for 
this increase is extracted from the thermal plasma, which leads to a 
small negative y-distortion with effective y-parameter. 


Up r ANy Op ANy 

yup(x.,Zi) ~ —Xi [l - e xiy/zi) —. (33) 


This counteracts the heating y-parameter, 


ctn r f AA'v 

yb(x„zd * [l - (34) 

caused by the BR absorption process. To fully include the effect of 
BR absorption, we simply need to multiply the scattering solution 
and yup by the survival probability given in Eq. (32). For x < 1, we 
thus have the Green’s function 


Gin(v, v',z) 


cpyiTo) e / [ln(x/x')-y^(z)]* 

-^— exp- 

47r .sjAnyy(z)x' \ 4y^(z) 


(1- 




> A 


(35) 


We find this approximation to work very well as long as corrections 
to the absorption optical depth caused by Doppler broadening are 
small (see Fig. 5). In particular, for Xi - 0.1 - 1 the solution works 
extremely well even until zi - 3 x 10*. 

The solution in Eq. (35) shows that, like in the yu-era, if pho¬ 
tons are injected only at very low frequencies, a high-frequency 
y-distortion appears through the net competition of heating (by BR 
absorption) and cooling (by low-frequency photon up-scattering). 
While at sufficiently low frequencies BR absorption can extract al¬ 
most all the injected photon energy, the cooling caused by scat¬ 
tering is limited to a small fraction oc y^ « 1. The transition fre¬ 
quency separating the regions of net heating to net cooling can be 



Figure 6. Comparison of the approximation in Eq. (38) with the full nu¬ 
merical results for xy = 5 and several injection redshifts. We also show the 
classical solution, Eq. (36), for z; = 3 X 10*, which clearly demonstrates the 
improvement of the new approximation. 


estimated with the condition 2yy{z) ~ T(f(Xh,z), as long as yy(z) is 
not too large. For 10* ;£ z < 10*, we find Xh - 0.01 - 0.1 (see 
Fig. 10), in very good agreement with our numerical calculations. 
At 0.01 < Xi < 1, the y-type contribution to the distortion caused 
by energy exchange and absorption remains relatively small. 


2.3.3 Photon injection at 10* < z ;£ 10* and 1 < x < 30 


To describe the solution at higher frequencies (1 < x < 30), we 
generally need to resort to numerical solutions. Neglecting recoil, 
one can use the classical solution (Zeldovich & Sunyaev 1969) 


An{x,y) = 


exp (-[ln(x/xi) - 3y]*/4y) 
yjAtry y’ 


(36) 


which differs from the low-frequency solution, Eq. (21), only by 
the net drift term (Chluba & Sunyaev 2008). The solution for pure 
recoil (neglecting any line broadening through recoil) simply is 
An(x,y) = A x“* 5[x-Xi(y)], with Xi(y) = Xi/(1-HXiy), which gives a 
drift Av/v == -xiy towards lower frequencies. One simple improved 
approximation, valid for Xjy « 1, thus is 


^ exp (- [ln(x/Xi) - 3y -I- ln(l -I- x,y)]* /4y) 

A„*(x,y)=^=-^^-*. (37) 

-jAny ^ 

This solution gives Apy(yy)lpy - (Apylpy)e*^y 1(1 + X[yy), which 
captures the aforementioned effects. 

We compared the numerical solution from simple diffusion 
calculations with this approximation and found that for larger val¬ 
ues of y and x-,, the position of the line was too low and the 
width a bit too large. Replacing the dispersion of the Gaussian 
by y — » y/(l -I- X\y) reproduced the width extremely well, even 
for larger values of y and X;. The match in the position of the 
line was further improved by replacing -3y —> -3y/ -^1 -l- xiy. 
To improve the match for Xi =: 1, we need to transition from 
-3y —» -y around x\ — 1. After several attempts, we found 
-3y —> -y[3 - 2/(xi)] with f(x\) = e^'^l + 2cf/2) to work very 
well. The match for the dispersion of the line was further improved 
by replacing y -» y/[l -l- x,y(l - /(x,))]. 
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Figure 7. Photon injection Green’s function for injection at intermediate redshifts, 5 x 10"^ Zi 3 x 10^. The photon injection Green’s function shows a rich 
phenomenology. To convert between x and observing frequency today we may use x ^ 0.017 (v/GHz). 
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Photon injection 11 


These considerations lead to the refined scattering solution 

exp (- [ln(x/Xi) - ay + ln(l + xj)]^ t^yP] 

An(x,y)=A -^- ;- -(38) 

yjAnypx^ 

with ft = [3 - 2f(x[)]/ and p - (I + Xiv[l - /(xi)])“*. 

The average energy density of the injected photons thus scales as 
Apy(yy)lpy = (Apy jpy) 1(1 + X'^yy). Followlng similar argu¬ 

ments as above, for x > 1 we then find 


Gm(v, V',Z) 


^-[ln{xlx')~ayy+\n{\+x'yy)fpyyP 

47r xjAnyjpx' 


-I- 



g4y.y(Q'+/3)g-Tff(y,s) 

1 + 



(39) 


where a and p are evaluated at x' and yy(z)- A comparison with 
the numerical results for xi = 5 and several injection redshifts is 
shown in Fig. 6. Clearly, the new approximation represents the full 
numerical results very well. 

Over a wider range of injection energies, Eq. (39) works very 
well until Zj 3 x 10^ (see Fig. 5). For x; 1 - 5 we found this 
solution to work even better, reaching up to Zi - 5 x 10“*. At high 
frequencies, photon absorption is already negligible and we can see 
from Fig. 5 that the net heating/cooling, which gives rise to a y-type 
contribution, can usually be neglected unless we inject at x; > 1 /Vy, 
for which recoil becomes significant. 

AtZi <: 3x10^*, the evolved line (omitting the y-part) no longer 
is well approximated by a simple Gaussian, with third moments 
becoming important (see Fig. 5). Improved approximations that in¬ 
clude higher order moments and frequency-dependent dispersion 
terms may be possible, but we leave this question to future work. 
In addition, closer to Zi - 5x10“*, corrections to the y-type contribu¬ 
tion due to the r-type (non-;r/non-y) distortion become significant. 
This could be captured by computing the effective heating rate from 
the evolution of the line as a function of redshift and then feeding it 
into the thermalization Green’s function of energy release to threat 
the heating contribution more precisely. 


2.4 The p-y transition era 

The signatures of photon injection during the p-y transition era 
(10‘*<z^3xl0^) show the richest phenomenology. In this regime, 
direct information about the initial distribution of photons can in 
principle be regained, since Comptonization is no longer able to 
smear photons out over the whole CMB energy spectrum like dur¬ 
ing the p-era. This is also the regime where heating of the matter 
by the injected photons becomes incomplete, so that the distortion 
starts to be dominated by the evolution of the injected photons when 
approaching zi —> 10“* and later (Sect. 2.3). 

In Fig. 5 and 7, we illustrate the numerical results for several 
cases. In particular, for injection at high frequencies (x; > 1) and 
Zi <: 5 X 10“*, the distortion shows large similarities with the dis¬ 
tortions from pure energy release. However, due to the addition of 
photons, the Green’s function for photon injection has a significant 
contribution oc G(v), especially when Zi —> 3 x 10^. This is be¬ 
cause the injected photons are smeared out over the whole CMB 
frequency range via Compton scattering without being strongly at¬ 
tenuated by photon absorption and a photon survival probability 
close to unity (see Fig. 8). We can also see that for x; = 1, a nega¬ 
tive y-type contribution arises because on average the plasma cools 
while smearing the injected photons out over the CMB spectrum. 
We find the transition between net heating and net cooling to occur 
at Xj 3.6 - 3.83, depending on the injection redshift (see Fig. 10). 



Figure 8. Survival probability for different injection frequencies and red- 
shifts after the /j-era. The curves were computed using CosmoTherm. At 
low redshifts (zi < 5 x 10“*), Compton scattering becomes inefficient, so that 
the survival probability steepens from 'Pa k er’‘''A m k ^ where xj 

can be deduced from Eq. (32). 


For injection at lower frequencies, in Fig. 7 one can still di¬ 
rectly identify the broadened and partially up-scattered photon line 
until the y-parameter exceeds unity significantly (zi > 2x 10^). This 
is because low-frequency photons have to Comptonize significantly 
until reaching the maximum of the CMB spectrum, a process that 
requires many scatterings. Focusing on the high-frequency distor¬ 
tion, for Xi = 10“^, one can also see the transition from net heating 
to net cooling, which occurs around zi - 2 x 10^ (see Fig. 10). 

In summary, the signals created by photon injection show a 
richer phenomenology than those caused by single energy release, 
in particular at zi £ 3 x 10^, where the final spectrum is found in a 
partially Comptonized state. If photons are injected at several fre¬ 
quencies, a superposition of dilferent distortion shapes can leave 
even richer signatures in the CMB spectrum. However, this also 
makes it harder to interpret the constraints on individual scenar¬ 
ios in a model-independent way, a problem that will be considered 
more carefully in a subsequent paper. 

3 PHOTON INJECTION AT HIGH ENERGIES 

The discussion of the preceding sections was limited to photon in¬ 
jection at X; < 30. Here, we consider injection at higher energies 
and Zi ^ 10^. For energies below the pair creation threshold with a 
soft background photon, Xp Im^PlkTy 4.3x 10*/(1-fz), the in¬ 
jected photons mainly transfer their energy to the medium via elec¬ 
tron recoil. In the expanding Universe, we have Xi(y) = Xi/(1 -l-xiy), 
so that the injected photon energy density is roughly given by 

Apyiy) a„Xi ANy 

-AAPPa ^ -r. (40) 

Py 1-1- Xjy Ny 

This approximation neglects any line broadening though electron 
recoil and Doppler terms, which are discussed in Sazonov & Sun- 
yaev (2000), but this should only lead to a correction. Equation (40) 
implies an energy release history 

d( 6 /Pr) „ cr-iNeC 

dz (l-l-Xiy)2 Ny H(\ + zY 
which can be directly used with the energy release Green’s func¬ 
tion to compute the distortion signal. Depending on the injection 
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epoch, eventually most of this injected photon energy heats the 
plasma, causing a /j-, y- and r-type distortion, with total energy re¬ 
lease ti.pylpy ^ O'pXi ANyINy. At eudy times (zi > fewx 10^, where 
Compton scattering is efficient) or for injection at 1/y^ x, < Xp 
(so that electron recoil is strong), the energy release for a single 
photon injection occurs quasi-instantaneously and can be estimated 
by adding Apylpy ^ UpX, ANyjNy at zi to the plasma. However, at 
lower redshifts, the shape of the distortion depends more directly 
on the energy loss rate. In this case, differences between the ‘on- 
the-spot’ approximation® and the detailed energy release history 
may be significant, especially when considering electromagnetic 
cascades above the pair-production threshold. 


3.1 Contributions from primary and secondary photons 


The picture given above is rough. In particular, it assumes that the 
contributions from secondary photons and non-thermal electrons, 
produced while the hard photons lose their energy, are negligi¬ 
ble. Similarly, the injected primary photon spectrum is neglected. 
This last statement is valid simply because for Xi » 1 and given 
Apylpy, below the COBE/FIRAS limit, Apyjpy < 6 x 10“® (Fixsen 
et al. 1996; Fixsen 2009), the added primary photon number is 
ANyjNy = iXp' xj^Apyjpy Apyjpy. Thus, ouly when N X\ sec¬ 
ondary photons are produced, does one have to worry about more 
than just the heating caused by the injected photon. 

One way for the primary photon to produce secondary photons 
is through DC emission with thermal electrons. Some of these DC 
photons are directly injected close to the CMB bands, x < 30. To 
estimate by how much this could change the final spectral distor¬ 
tion, we shall use the DC emissivity in the soft photon limif (Light- 
man 1981; Thorne 1981); however, we add relativistic corrections 
due to the fact that the injected photon can be hard, hvi « m^c^. 
This gives (Chluba et al. 2007b) 
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We neglected stimulated DC emission due to the scattered photon 
y' Vi, which should cause a tiny correction. We furthermore de¬ 
fined the temperature ratio, (j>c(y) = TyjTciy), where Tdy) is the 
color temperature of the high-frequency photon, which initially ex¬ 
ceeds the CMB photon temperature by a large amount. 

Inserting n^^jx) = INy)x^^6[x — Xi(y)] for the distribu¬ 

tion of the hard photon, with Xj, = Xi/(1 +Xiy), we have the emission 
coefficient 
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For the CMB blackbody photons, most of the DC emission is pro¬ 
duced at X = (Jdc/G^y^^ - 3.3, with the effective DC integral 
-^dc = f x‘'«b(l -I- Mbb) ~ 25.976. Thus, the relative efficiency of DC 
emission caused by the hard photons is 
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(44) 


® This means that the injected photon energy is converted into heat quasi- 
instantaneously (Padmanabhan & Finkbeiner 2005). 

^ In this limit, the emitted DC photon has a much lower energy than the 
scattering high-energy photon. 


We find that even if we saturate the COBE/FIRAS bound, 
ANyjNy < 6 X 10“®/(apXi) « 1.6 x lO^^'xr*, for X[ < Xp, this extra 
emission relative to the CMB hlackbody DC emission can only be 
important at* Zi ^ 4000. Since at these low redshifts, DC is already 
much slower than BR, this only adds a small correction. This state¬ 
ment also holds when going beyond the soft photon limit (Gould 
1984; Chluba 2005), which allows including the DC emission at 
higher frequencies. 

How many secondary photons could maximally be produced? 
For x; below the pair-production threshold, the hard photons will 
indeed lose their energy mainly through recoil (Zdziarski & Svens- 
son 1989). At high redshifts (z > 3300), the Thomson optical depth, 
T « 0.21(1 -I- z), strongly exceeds unity, showing that photons un¬ 
dergo many scatterings. The loss of energy by recoil in a single 
scattering event (t 1) is Avjv - -hvjnicC^. Thus, as long as 
hVi <K TMeC^, the scattered thermal electron only receives a small 
kick, subsequently sharing its energy with the thermal plasma on a 
short time-scale. In this case, Eq. (41) describes the situation very 
well, since hardly any energy is converted into energetic secondary 
particles that could also enhance the BR and DC emissivities. 

For hVi > nicC^, a significant population of secondary, non- 
thermal electrons can be built up. These cool down by inverse 
Compton scattering with CMB photons, which in addition removes 
photons from the CMB bands, creating additional hard photons. 
This is equivalent to a photon destruction event (negative photon in¬ 
jection), which could leave a signature in the final distortion. How¬ 
ever, the total number of all energetic secondary particles (includ¬ 
ing up-scattered CMB photons) can never become dramatic due 
to energetic constraints. The only way to significantly enhance the 
number of photons may be by very soft photon production (x « 1), 
for example, through non-thermal BR. This does not require a lot 
of energy but could directly help in the thermalization process. 

To include these aspects consistently, more detailed cas¬ 
cade calculations in the pre-recombination era are required. These 
should not only address how the medium is heated as a function 
of time, an aspect that determines the exact distortion shape (p-, y- 
and r-distortion) caused by heating, but also how many secondary 
photons are eventually added to or removed from the CMB bands 
and when. At Zi £ few x 10®, the final secondary photon population 
reaching the CMB bands from the high energies may also leave its 
features. Combining these aspects could allow discerning between 
different scenarios related to decaying or annihilating particles. 


4 CURRENT CONSTRAINTS ON PHOTON INJECTION 

In this section, we give a brief discussion of current constraints 
on photon injection from COBE/FIRAS (Fixsen et al. 1996; Fixsen 
2009) and BBN (Simha & Steigman 2008; Jeong et al. 2014). From 
Simha & Steigman (2008), we have 

ANy 

^ -0.08 ± 0.07 (68% c.l.) (45) 

Ny 

a number that is derived by comparing the CMB temperature at re¬ 
combination and BBN, where the latter is derived using measured 
light element abundances, implying |AA,,/Ay| < 0.07 (68% c.l.). 
While this limit is not very tight, it supersedes the constraint from 
p < 9 X 10“® (95% c.l.) (Fixsen et al. 1996) for photon injection 


“ Neglecting the suppression of DC emission due to relativistic corrections, 
one obtains zi S 7 x 10**. 
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Figure 9. COBE/FIRAS limits (68% c.l.) on photon production for red- 
shifts during the /j-distortion era. Around JCi 3.6, the net chemical po¬ 
tential created by photon production vanishes so that the tightest limit is 
|AAy/Ay| < 0.07 (68% C.L) from measurements of the light element abun¬ 
dances (Simha & Steigman 2008). 


at very low frequencies and around x\ 3.6 (see Fig. 9). At high 
frequencies, one also has 


^Ny 

~n; 


< 


8.9 X 10^^ 

te-3.6)jr4 


(68% C.L), 


(46) 


which imposes a very tight limit on energetic photon production 
during the ^-era. 

To derive constraints on photon injection at different frequen¬ 
cies and for Zi £ 3 x 10^, one would need to consider the precise 
shape of the final distortion in light of COBE/FIRAS data, since 
the distortion can no longer be parametrized as a simple /r or y- 
distortion (see Fig. 5 and 7). In particular, for 1 < x; < 10, cor¬ 
responding roughly to the COBE/FIRAS channels 68 GHz < v < 
640 GHz, one could expect residual narrow features from the direct 
photon injection event at these redshifts. While we leave a detailed 
discussion for the future, for injection at very high and very low fre¬ 
quencies (outside the regime directly probed with COBE/FIRAS) 
one can still obtain some constraints. For injection at low frequen¬ 
cies (xi « 1), the overall heating of the plasma (which would leave 
a high-frequency yu/y-distortion) remains small unless we allow for 
a very large photon production, ANyINy == 1. In this case, more 
detailed computations of the distortion may be required, since non¬ 
linear effects could become important. This statement may also ap¬ 
ply to soft photon injection at high redshift. 

For injection of energetic photons (x; > 10) at zi < 3 x 10^, 
from Eq. (40) we have the total energy release 


Apy OpX^yy ANy 
Py 1 - 1 - X\yy Ny 


(47) 


between z\ and today. For photon injection at late times (10^ £ Zi £ 
5 X 10"^ and <K 1), this gives 


AW 


Ny 


8.1 X 10-5 

x^yy 


(68% C.L). 


(48) 


This limit weakens strongly towards lower redshifts, since only a 
small fraction oc xjy^ of the total injected photon energy is con¬ 
verted into heating of the medium, simply because electron recoil 
becomes inefficient. 
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Figure 10. Regions for net heating and cooling by photon injection. At high 
frequencies, depending on the injection redshift, the transition between net 
cooling and heating occurs around x; ^ 3.6 - 3.83, while at low frequencies 
it occurs over a wide range of frequencies at different redshifts (blue dash- 
dotted lines). In black, we also show the injection frequency for which the 
total recoil exceeds unity and the critical frequency at which Ps ~ 1/e 
(optical depth r « 1). The horizontal lines also indicate redshifts of constant 
Vy (dashed red) and constant distortion visibility (dashed purple). 

5 CONCLUSIONS 

We studied the spectral distortions signatures caused by photon in¬ 
jection over a wide range of energies and redshifts. Our calculations 
illustrate the rich phenomenology of the final distortion shapes, 
with several examples shown in Fig. 5 and 7. They also constitute 
the starting point for future investigations of distortions created by 
specific photon injection scenarios. 

In Sect. 2, we discussed a detailed picture for the shape of the 
photon injection Green’s function in different regimes. Figure 10 
summarizes the different domains. Generally, we find that for pho¬ 
ton injection at high (x; > 3.6 - 3.83) and very low frequencies 
(x; < IQ-'* - 10-5) ^ jjgj jjeating of the plasma occurs. At low fre¬ 
quencies this is because photons are efficiently absorbed and thus 
directly converted into heat, while at high frequencies it is because 
on average photons have to down-scatter to thermalize. At interme¬ 
diate regimes, a net cooling of the plasma occurs. 

For photon injection during the /r-era (z > 3 x 10^), simple an¬ 
alytic approximations for the amplitude of p are given in Sect. 2.2. 
Similarly, for moderate scattering y-parameter (zi £ 3 x 10“'), we 
found simple approximations for the photon injection Green’s func¬ 
tion. In particular, we obtained a new Compton scattering solution, 
Eq. (38) for single photon injection, which is valid over a wide 
range of energies. The solution merges the high- and low-frequency 
domain, simultaneously including Doppler broadening and boost¬ 
ing, electron recoil and stimulated scattering. At intermediate red¬ 
shifts, detailed numerical calculations are required, since the final 
distortion can be found in a partially Comptonized state (Fig. 5 
and 7). The introduced Green’s function method provides a simple 
way to accelerate the computations in this regime and will be made 
available at www. Chluba. de/CosmoTherm. 

Our calculations show that photon injection can offset the dis¬ 
tortion signals created by pure energy release. In particular, since 
photon injection can create negative p and y distortion contribu¬ 
tions, this may lead to significant limitations for the interpretation 
of future CMB spectral distortions measurements. However, for 

































14 Chluba 


photon injection at z; < 3 x 10^, partial direct information about the 
photon injection process may be recovered. This could help dis¬ 
tinguishing pure energy release and photon injection scenarios, a 
problem that will be considered in a future analysis. 

We briefly discussed photon injection at high energies, reach¬ 
ing close to and above the pair-production threshold (Sect. 3). 
Our analysis indicates that above the pair-production threshold, 
more detailed computations of the soft photon production efficiency 
(x < 1) may be required. These soft photons could help thermaliz- 
ing the distortion and also may introduce features into the final sig¬ 
nal that could further help discerning different processes. Similarly, 
detailed considerations of the reprocessing of photons and energy 
during the post-recombination by atomic and molecular species 
should be considered more carefully. 

In Sect. 4, we presented COBE/FIRAS and BBN constraints 
on photon injection. These could be improved significantly with 
a PIXIE-like experiment. The constraints in the yu-era are summa¬ 
rized in Fig. 9, showing that scenarios with significant photon pro¬ 
duction around Xj 0.01 - 10 are already excluded. However, at 
very low frequencies (Xj < 0.01) and around Xj 3.6, in principle 
significant photon injection (ANyjNy - 1) is still allowed and may 
require the inclusion of non-linear effects. For Xj < 0.01 this is sim¬ 
ply because the added energy Apyjpy ^ a^x^ANylNy is very small 
unless ANy/Ny becomes large, while for xy - 3.6 a balanced photon 
injection scenario, with vanishing net distortion (see Sect. 2.1.2), 
is encountered. For photon injection at lower redshifts, a detailed 
comparison with COBE/FIRAS or future spectral distortion data 
has to be carried out, since in principle narrow features can remain 
in the partially Comptonized regime. We leave this to future work. 

We close by noting that further strong motivation for consid¬ 
ering scenarios with photon injection can be given due to the ob¬ 
served low-frequency excess at hv 3.3 GHz (Fixsen et al. 2011), 
which may have an interpretation as a signature of a decaying or an¬ 
nihilating particle in the pre-recombination era at zi £ 10^. In this 
case, soft secondary photons could be produced that now are found 
in a partially Comptonized state, however, more detailed computa¬ 
tions for specific injection scenarios are required. 
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